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Abstract 

Ordinary approach to quantum algorithm is based on quantum Turing 
Oh, machine or quantum circuits. It is known that this approach is not pow- 

erful enough to solve NP-complete problems. In this paper we study a new 
£h ' approach to quantum algorithm which is a combination of the ordinary 

quantum algorithm with a chaotic dynamical system. We consider the 
satisfiability problem as an example of NP-complete problems and argue 
that the problem, in principle, can be solved in polynomial time by using 
our new quantum algorithm. 
keywords: Quantum Algorithm, NP-complete problem, Chaotic Dynam- 



1 Introduction 

Ordinary approach to quantum algorithm is based on quantum Turing machine 
or quantum circuits 0121 EH ■ It is known that this approach is not powerful 
enough to solve NP-complete problems 0J |S] . In we have proposed a new 
approach to quantum algorithm which goes beyond the standard quantum com- 
putation paradigm. This new approach is a sort of combination of the ordinary 
quantum algorithm and a chaotic dynamics. This approach was based on the 
results obtained in the paper [7|. 

There are important problems such as the knapsack problem, the travel- 
ing salesman problem, the integer programming problem, the subgraph isomor- 
phism problem, the satisfiability problem that have been studied for decades 



and for which all known algorithms have a running time that is exponential in 
the length of the input. These five problems and many other problems belong 
to the set of NP-complete problems 0]. 

Many NP-complete problems have been identified, and it seems that such 
problems are very difficult and probably exponential. If so, solutions are still 
needed, and in this paper we consider an approach to these problems based on 
quantum computers and chaotic dynamics as mentioned above. 

As in the previous papers [71IHJ we again consider the satisfiability problem as 
an example of NP-complete problems and argue that the problem, in principle, 
can be solved in polynomial time by using our new quantum algorithm. 

It is widely believed that quantum computers are more efficient than classical 
computers. In particular Shor |H||U] gave a remarkable quantum polynomial-time 
algorithm for the factoring problem. However, it is known that this problem is 
not NP-complete but is NP-intermidiate. 

Since the quantum algorithm of the satisfiability problem (SAT for short) 
has been considered in jJJ, Accardi and Sabbadini showed that this algorithm 
is combinatric one and they discussed its combinatric representation |10|. It 
was shown in [7] that the SAT problem can be solved in polynomial time by 
using a quantum computer under the assumption that a special superposition 
of two orthogonal vectors can be physically detected . The problem one has to 
overcome here is that the output of computations could be a very small number 
and one needs to amplify it to a reasonable large quantity. 

In this paper we construct a new model (representation) of computations 
which combine ordinary quantum algorithm with a chaotic dynamical system 
and prove that one can solve the SAT problem in polynomial time. 

For a recent discussion of computational complexity in quantum computing 
see ^2 El E3 EU Mathematical features of quantum computing and quantum 
information theory are summarized in |15| . 

2 SAT Problem 

Let X = {xi, • • • , x n } be a set. Then Xk and its negation Xk (k = 1, 2, • • • , n) are 
called literals and the set of all such literals is denoted by X — {x\,x\, • • • , x n ,x n }. 
The set of all subsets of X is denoted by T{X ) and an element C € T[X ) is 
called a clause. We take a truth assignment to all variables Xk- If we can assign 
the truth value to at least one element of C, then C is called satisfiable. When 
C is satisfiable, the truth value t(C) of C is regarded as true, otherwise, that of 
C is false. Take the truth values as true "1", false "0". Then 

C is satisfiable iff t(C) = 1. 

Let L — {0, 1} be a Boolean lattice with usual join V and meet A, and t(x) 
be the truth value of a literal x in X. Then the truth value of a clause C is 
written as 

t(C) = \Z xe ct(x). 



Further the set C of all clauses Cj (j — 1, 2, • • • m) is called satisfiable iff the 
meet of all truth values of Cj is 1; 

t(C) = A^itfC,-) = 1. 

Thus the SAT problem is written as follows: 

Definition 1 SAT Problem: Given a set X = {x\, • • • , x n } and a setC = {C\, Ci, • • • , C m } 
of clauses, determine whether C is satisfiable or not. 

That is, this problem is to ask whether there exsits a truth assignment to 
make C satisfiable. 

It is known0] in usual algorithm that it is polynomial time to check the satis- 
fiability only when a specific truth assignment is given, but we can not determine 
the satisfiability in polynomial time when an assignment is not specified. 

Note that a formula made by the product (AND A) of the disjunction (OR 
V) of literals is said to be in the product of sums (POS) form. For example, the 
formula 

(x 1 Vi 2 ) A (xi) A (i 2 Vi 3 ) 

is in POS form. Thus a formula in POS form is said to be satisfiable if there is 
an assignment of values to variables so that the formula has value 1. Therefore 
the SAT problem can be regarded as determining whether or not a formula in 
POS form is satisfiable. 

The following analytical formulation of SAT problem is useful. We define a 
family of Boolean polynomials fj., indexed by the following data. One A is a 
set 

A = {Si, ..., Sn,Ti, ...,T/v} , 

where Si,Ti C { 1 , . . . , n} , and /^ is defined as 

j/U (Si,*" ,X n ) ■- 




i=l 



We assume here the addition modulo 2. The SAT problem now is to determine 
whether or not there exists a value of x = (x\, ■ ■ ■ , x n ) such that /a(x) = 1. 

3 Quantum Algorithm 

Although the quantum algorithm of SAT problem is needed to add the dust 
bits to the input n bits, the number of dust bites has been shown the order of n 
[7|E3. Therefore for simplicity we will work in the (n + l)-tuple tensor product 
Hilbert space TL = ®" +1 C 2 in this paper with the computational basis 

\xi,...,x n ,y) =®i =1 \Xi)<8) \y) 



where Xi,...,x n , y = or 1. We denote \x\, ...,x n ,y) = |x, y) . The quan- 
tum version of the function /(x) := /a(x) is given by the unitary operator 
Uf |x, y) = |x, y + /(x)) . We assume that the unitary matrix Uf can be build 
in the polynomial time, see [7|- Now let us use the usual quantum algorithm: 
(i) By using the Fourier transform produce from |O,0) the superposition 
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(ii) Use the unitary matrix Uf to calculate /(x) : 

* X 

Now if we measure the last qubit, i.e., apply the projector P = I®|l)(l|to 
the state \vf) , then we obtain that the probability to find the result /(x) = 1 
is \\P \vf)\\ — r/2 n where r is the number of roots of the equation /(x) = 1. If 
r is suitably large to detect it, then the SAT problem is solved in polynominal 
time. However, for small r, the probability is very small and this means we in 
fact don't get an information about the existence of the solution of the equation 
/(x) = 1, so that in such a case we need further deliberation. 

Let us simplify our notations. After the step (ii) the quantum computer will 
be in the state 

\v f ) = v 7 !-? 2 \<Po) <8> |0) + q \(pi) ® |1) 



where \(p\) and \(fo) are normalized n qubit states and q = yr/2". Effectively 
our problem is reduced to the following 1 qubit problem. We have the state 



and we want to distinguish between the cases q — and q > (small positive 
number). 

It is argued in [5] that quantum computer can speed up NP problems 
quadratically but not exponentially. The no-go theorem states that if the inner 
product of two quantum states is close to 1, then the probability that a mea- 
surement distinguishes which one of the two it is exponentially small. And one 
could claim that amplification of this distinguishability is not possible. 

At this point we emphasize that we do not propose to make a measurement 
(not read) which will be overwhelmingly likely to fail. What we do it is a 
proposal to use the output \ip) of the quantum computer as an input for another 
device which uses chaotic dynamics in the sequel. 

The amplification would be not possible if we use the standard model of 
quantum computations with a unitary evolution. However the idea of our paper 
is different. We propose to combine quantum computer with a chaotic dynamics 
amplifier. Such a quantum chaos computer is a new model of computations 
going beyond usual scheme of quantum computation and we demonstrate that 
the amplification is possible in the polynomial time. 



One could object that we don't suggest a practical realization of the new 
model of computations. But at the moment nobody knows of how to make 
a practically useful implementation of the standard model of quantum com- 
puting ever. Quantum circuit or quantum Turing machine is a mathematical 
model though convincing one. It seems to us that the quantum chaos computer 
considered in this paper deserves an investigation and has a potential to be 
realizable. 

In this paper we propose a mathematical model of computations for solving 
SAT problem by refining our previous paper |Hj. A possible specific physical 
implementation of quantum chaos computations with some error correction will 
be discussed in a separate paper |16| . which is some how related to the recently 
proposed atomic quantum computer |17| . 

4 Chaotic Dynamics 

Various aspects of classical and quantum chaos have been the subject of numc- 
rious studies, see JH1 and ref's therein. The investigation of quantum chaos by 
using quantum computers has been proposed in ^3 [23 IS] • Here we will argue 
that chaos can play a constructive role in computations. 

Chaotic behaviour in a classical system usually is considered as an exponen- 
tial sensitivity to initial conditions. It is this sensitivity we would like to use to 
distinquish between the cases q — and q > from the previous section. 

Consider the so called logistic map which is given by the equation 

X n +l = ax n (l - X n ) = g{x), X n € [0, 1] . 

The properties of the map depend on the parameter a. If we take, for example, 
a = 3.71, then the Lyapunov exponent is positive, the trajectory is very sensitive 
to the initial value and one has the chaotic behaviour |18) . It is important to 
notice that if the initial value xo — 0, then x n = for all n. 

It is known |5j that any classical algorithm can be implemented on quantum 
computer. Our quantum chaos computer will be consisting from two blocks. 
One block is the ordinary quantum computer performing computations with the 
output \ip) = \J\ — q 2 |0) + q\l). The second block is a computer performing 
computations of the classical logistic map. This two blocks should be connected 
in such a way that the state \i/j) first be transformed into the density matrix of 
the form 

p = q 2 Pi + (1 - q 2 ) P 

where P\ and Pq are projectors to the state vectors |1) and |0) . This connection 
is in fact nontrivial and actually it should be considered as the third block. 
One has to notice that P\ and Pq generate an Abelian algebra which can be 
considered as a classical system. In the second block the density matrix p above 
is interpreted as the initial data po, and we apply the logistic map as 

_ _(I + g m ( P o)a 3 ) 



where I is the identity matrix and 173 is the z-component of Pauli matrix on C 2 . 
This expression is different from that of our first paper 0. To find a proper 
value m we finally measure the value of 173 in the state p m such that 

M m = trp m a 3 . 
After simple computation we obtain 

Pm= (I + 9^) ^ dMm = gm{q2) 

Thus the question is whether we can find such a m in polynomial steps of 
n satisfying the inequality M m ^ i for very small but non-zero q 2 . Here we have 
to remark that if one has q = then p — P Q and we obtain M m = for all m. If 
q ^ 0, the stochastic dynamics leads to the amplification of the small magnitude 
q in such a way that it can be detected as is explained below. The transition 
from po to p m is nonlinear and can be considered as a classical evolution because 
our algebra generated by Pq and Pi is abelian.The amplification can be done 
within atmost 2n steps due to the following propositions. Since g m (q 2 ) is x m 
of the logistic map x m +i — g(x m ) with xo — q 2 , we use the notation x m in the 
logistic map for simplicity. 

Proposition 2 For the logistic map x n+ \ — ax n (1 — x n ) with a £ [0,4] and 
xq £ [0, 1] , let xq be ^r and a set J be {0, 1, 2, • • • , n, • • • 2n} . If a is 3.71, then 
there exists an integer m in J satisfying x m > i . 

Proof: Suppose that there does not exist such rn in J. Then x m < \ for any 
m £ J. The inequality x m < \ implies 



6.lL[l Xm—ljXm — l r_ X m — i. 



Thus we have 



1 3.71 /3.71Y /3.71Y 1 

2- Xm ~— Xm - 1 -■■■-{—) Xa ={—) ^ 
from which we get 

2 n+m-i > (3 71 ) m _ 

According to the above inequality, we obtain 

n-1 

m < 



log 2 3.71 - 1 ' 
Since log 2 3.71 4= 1.8912, we have 

n — 1 5 . . 

m < < — (n — 1 , 

~ log 2 3.71 - 1 4 v ' 

which is definitely less than In — 1 and it is contradictory to the statement 
"x m < i for any m £ J". Thus there exists m in J satisfying x m > i.B 



Proposition 3 Let a and n be the same in the above proposition. If there exists 



1 +h„r, ™_ ^ »»— l 



rrio in J such that x mo > ■* , then too > -^ — wri' 
Proof: Since < x m < 1, we have 

•Em = o. t lyi X m — i)X m —i \ o. t LXm—Xi 

which reduces to 

x m < (3.71) m x . 

For too in J satisfying x mo > ^ , it holds 

1 1 



(3.71) m °'"° 2x (3.71) r 



It follows that from xo 



which implies 



l 



log 2 2 x (3.71) mo > n 



n — 1 
too > 



log 2 3.71 

According to these propositions, it is enough to check the value x m (M m ) 
around the above too when q is -^ for a large n. More generally, when q=^ 
with some integer k, it is easily checked that the above two propositions are 
held and the value x m (M m ) becomes over ^ around the Too above. 

One can think about various possible implementations of the idea of using 
chaotic dynamics for computations, which is open and very intersting problem. 
About this problem, realization of nonlinear quantum gates will be essential, on 
which we will discuss in atomic quantum computer in |16|. 

Finally we show in Fig.l how we can easily amplify the small q in several 
steps. 



5 Conclusion 

The complexity of the quantum algprithm for the SAT problem has been con- 
sidered in [x] where it was shown that one can build the unitary matrix Uf in 
the polynomial time. We have also to consider the number of steps to in the 
classical algorithm for the logistic map performed on quantum computer. It is 
the probabilistic part of the construction and one has to compute several times 
to be able to distingish the cases q = and q > 0. Thus it concludes that 
the quantum chaos algorithm can solve the SAT problem in polynominal time 
according to the above propositions. 

In conclusion, in this paper the quantum chaos algorithm is proposed. It 
combines the ordinary quantum algorithm with quantum chaotic dynamics am- 
plifier. We argued that such a algorithm can be powerful enough to solve the 



NP-complete problems in the polynomial time. Our proposal is to show exis- 
tence of algoritm to solve NP-complete problem. The physical implimentation 
of this algorithm is another question and it will be strongly desirable to be 
studied. 
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